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throughout  the  axial  distance  covered  by  these  tests.  The  jet-prcpelled 
body  produces  the  middle  sized  wake  compared  to  the  other  two  models  up  to 
station  X/D  = 20;  beyond  this  station,  this  wake  becomes  the  smallest.  In 
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the  highest  Reynolds  shear  stresses  near  the  tail  of  the  body.  Because 
this  wake  exhibits  the  fastest  decay  rate  for  all  the  turbulence  properties 
measured  it  has  the  lowest  values  of  these  properties  at  the  farthest 
station  (X/D  = 40). 

The  results  are  compared  with  the  available  data  obtained  by  other 
researchers.  Momentum  and  energy  analyses  are  used  to  check  the  accuracy 
of  the  results  and  to  obtain  an  insight  into  the  flow  structure  and  the 
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I.  INTRODUCTION 


Tha  turbulent  wake  behind  a body  in  a fluid  stream  has  been  a 
frequent  subject  of  experimental  research  because  the  development  of 
such  flow  is  not  fully  understood,  and  the  turbulent  properties 
(Reynolds  shear  stress  and  turbulence  intensity)  can  only  be  calculated 
with  the  aid  of  semiempirical  assumptions.  Several  experimental 
studies  have  been  made  at  low  and  high  speeds  for  both  axisymmetric  and 
planar  geometries.  Table  I provides  a summary  of  the  more  imoortant 
work  concerned  with  subsonic,  turbulent  wakes.  Schlichting  [1], 

Page  and  Falker  [2],  Hall  and  Hislop  [3],  and  Hill  et  al  [4]  measured 
the  mean  flow  properties  in  the  wake  behind  cylinders.  Townsend  [5] 
obtained  some  data  on  the  turbulent  properties  of  the  cylinder  wake. 
Cooper  and  Lutzkv  [6],  and  Carmodv  [7]  used  their  disc  models  and  mea- 
sured both  the  mean  flow  and  turbulent  characteristics.  More  recently, 
wake  studies  have  been  made  for  an  ellipsoid  18],  a square  cylitder  [9], 
a spheroid  [10],  and  a streamlined  body  [11]. 

During  the  past  few  years,  turbulent  wake  flows  with  near-zero 
momentum  defect  have  become  a topic  of  increasing  interest  to 
researchers.  Such  flows  occur  behind  many  self-propelled  vehicles  and 
devices,  including  submarines,  rockets,  missiles,  etc.  An  under- 
standing of  near-zero  lomentum  (or  momentumless)  wake  flows  is  neces- 
sary in  order  to  approach  several  practical  problems  such  as  submarine 
detection,  |>oliut<on  dispersion  from  propulsion  units,  and  wake  dis 
si  pa  I ion.  Only  verv  tew  works  .haling  with  the  .iv  ment  un.  I . . war... 

soil -propel  led  bodies  have  been  published,  namely  those  by 
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Summary  of  Detailed  Measurements  of  the  Vain  ''ixing  Region  in  Turbulent  V’akes 
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Naudascher [12]  Ginevskii  et  al  [13], and  Gran[l4].  The  configuration  used 

by  Naudascher  and  Ginevskii  was  a circular  disk  normal  to  an  oncoming 

fluid  stream.  A higl  speed  axial  jet  exhausted  from  the  disc  center 

was  applied  to  enforce  a zero  net  momentum  wake.  The  wind  tunnel 

velocity  and  Reynolds  number  in  Naudascher 's  tests  were  60  ft/sec,  and 

4 

5.4x10  , respectively,  and  single  and  cross  hot  wives  were  used  to 
measure  the  turbulence  characteristics.  The  models  used  by  Naudascher 
and  Gir.evskii  could  not  give  sufficiently  good  simulation  of  most  of 
the  cases  of  practical  interest  becav^e  their  models  were  neither 
slender  nor  streamlined  as  is  the  usual  case  of  self-propelled  bodies. 
Gran  studied  the  momentumless  wake  of  a slender,  propeller-driven  body 
having  a Rsnkine  ovoid  contour  (Ref . 14) . His  measurements  included 
both  the  mean  flow  and  turbulent  properties  of  the  wake.  The  Reynolds 
number  was  relatively  low  (ReD  ~ 6 x 10*)  and  some  difficulties  with 
obtaining  symmetrical  wakelwere  encountered.  As  Table  I demonstrates, 
most  of  the  available  information  is  related  to  unstreamlined  bodies 
and  was  obtained  under  the  conditions  of  low  Reynolds  numbers.  More- 
over, only  two  distinctly  different  moraentumless  cases  have  been 
investigated  thus  far. 

In  view  of  the  scarcity  and  limitations  of  the  available  informa- 
tion on  turbulent  momentumless  wakes  it  was  felt  that  an  experimental 
investigation  was  in  order  that  would  cover  more  practical,  realistic 
wakes  and  possibly  high  Reynolds  number  range.  The  program  was  con- 
ceived to  provide  a systematic  comparison  of  the  turbulent  wakes 
behind  three  separate,  but  related,  slender  bodies:  (1)  a streamlined 
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drag  body  (Model  I),  (2)  jet-propelled  body  (Model  II),  and  (3) 

propeller-driven  body  (Model  III).  All  bodies  bad  identical  blunt 
noses  and  parallel  middle-bodies  and  all  had  streamlined,  sharp  stern 
sections.  The  maximum  diameter  of  each  of  the  three  models  was  six 
inches  and  the  length-to-diameter  ratio  was  12  for  Models  I and  III 
and  10.29  for  Model  II.  The  experiments  were  carried  out  in  the 
VPI&SU  Stability  Wind  Tunnel  at  a speed  of  206  ft/sec.  The  Reynolds 
number  based  on  diameter  was  6.18  x 10"*,  i.e,  an  order  of  magnitude 
larger  than  for  previous  propeller-driver  laboratory  tests.  For  each 
model,  the  wake  properties  were  measured  at  five  stations  (X/D  =2,  5, 
10,  20,  and  40).  Mean  flow  properties  were  determined  with  a pitot- 
static  probe,  a wedge  type  directional  probe,  and  a Kiel  probe.  The 
following  turbulence  intensity  and  Reynolds  shear  stress  values  were 
obtained  using  single-normal  and  single-yawed  hot-wire  probes: Vu^2, 
Vv'2,  Vw'2,  u'v',  and  u'w' . Momentum  and  energy  balance  analyses 
were  carried  out  to  verify  the  measurement  accuracies  and  to  obtain 
some  insight  into  the  mechanisms  of  momentum  and  energy  transformation. 
In  this  work  only  the  turbulence  date  are  evaluated  in  detail  and  com- 
pared with  the  experimental  results  of  Chevr^y,  Waudascher  and  Gran. 

The  mean  flow  data  are  presented  and  discussed  in  Ref.  19. 
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II.  WIND  TUNNEL  AND  MODELS 


Wind  Tunnel 

The  experiments  were  performed  in  the  VPI  stability  wind  tur.nel 
(Fig.  1)  at  a free  stream  velocity  corresponding  to  9.5  in.  H^O  dyna- 
mic pressure  (about  206  ft/sec).  The  facility  has  a test  section 
6 ft  x 6 ft,  23  ft  long,  and  the  cross  section  diverges  slightly  in 
the  downstream  direction.  The  mean  velocity  at  the  edge  of  the  wake 
for  the  five  stations  varied  less  than  6%,  and  the  free  stream  tur- 
bulence intensity  outside  the  wake  was  measured  to  be  less  than  0.1Z. 

Models 

Three  separate,  but  related,  models  were  tested: 

(i)  Model  I - Streamlined  Drag  Body  (Fig.  2) 

Model  I is  a body  of  continuous  slope  variation  with  a parabolic 
nose  made  of  wood,  a 3 ft  long  section  of  parallel  middle  body  made  of 
aluminum  tube,  and  a sharp  wood  tail.  The  total  length  of  the  body  is 
72  inches. 

(ii)  Model  II  - Jet  Propelled  Body  (Figs. 3 and  4) 

Model  II  is  propelled  by  a peripheral  air-jet.  It  has  a parabolic, 
plexiglass  nose,  a parallel  aluminum  middle  section,  and  a stern  which 
is  formed  of  a plexiglass,  tear  drop  shaped  center  body  supported  by 
four  brass  vanes  spaced  90  degrees  apart.  Air  was  injected  through  a 
1 inch  peripheral  slot  at  the  end  of  the  parallel  middle  section.  The 
injection  was  supplied  by  four  150  psi  compressors;  a control 
valve  was  used  to  adjust  the  pressure  of  the  injected  air  (about 
112  psia)  so  that  the  wake  was  approximately  momentumless. 
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The  air  was  fed  by  a 3 in.  diameter  pipe  outside  of  the  test  section 
and  a supporting  strut  (made  of  seventeen  0.5-in.  diameter  pipes)  into 
a planum  chamber  inside  the  model.  The  plenum  chamber  was  vented 
through  a series  of  coarse  and  fine  screens  a:.d  exhausted  through  the 
slot.  The  screens  were  used  to  minimize  nonuniformities  and  turbulence 
of  the  exhausted  jet. 

(iii)  Model  III  - Propeller-Driven  Body  (Fig.  5) 

The  propeller-driven  model  was  actually  the  drag  boay  model  modi- 
fied so  that  a propeller  and  an  electric  motor  could  be  installed.  A 
model  2M145  Dayton,  1.2  hp , 10,000  rpm,  AC/DC  motor  was  fitted  inside 
the  model,  and  a 6-in.  diameter  propelle-  was  connected  to  the  motor 
using  a 23-3/8  inch  long  shaft.  The  rpm  was  measured  by  a magnetic 
pick-up  device.  The  signal  produced  by  the  magnetic  pick-uo  was  first 
amplified  and  then  displayed  on  the  oscilloscope  screen.  The  propeller 
was  a three-bladed , plastic  model  airplane  propeller.  In  order  to 
increase  the  pitch  angle  of  the  blades,  they  were  heated  and  twisted 
by  using  a special  tool  made  for  this  purpose.  The  experimental  mean 
diameter  pitch  of  the  propeller  was  about  2.5,  and  the  propeller  was 
operated  at  an  advance  ratio  J = 2 with  an  apparent  efficiency  68%. 

The  de' ailed  performance  of  the  propeller  was  reported  by  Swanson  [19]. 

To  obtain  a raomentumless  wake,  the  propeller  had  to  be  driven  at 
speeds  between  12,000  and  12,500  rpm,  i.e.,  at  speeds  higher  than  the  motor 
was  designed  for.  This  resulted  in  some  overloading  of  the  motor.  In 
order  to  prevent  excessive  temperature  rise,  the  unit  required 
an  efficient  water  cooling  system.  The  latter  was  provided  by  a 
svstem  of  coils  made  of  1/2  in.  copper  tubing  and  brazed  to  the  motor 
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case.  A thermocouple  was  hooked  up  to  monitor  the  temperature  of  the 
motor,  and  this  served  to  check  whether  the  motor  was  functioning 
properly. 

Model  I anj  III  were  suspended  from  the  ceiling  of  the  wind 
tunnel  by  a thin  strut,  made  of  aCl5  inch  thick  steel  plate.  There 
were  0.5  in.  diameter  pipes  attached  to  the  leading  edge  of  the  plate 
and  they  served  as  conduits  for  the  electrical  wiring  and/or  as  chan- 
nels for  the  cooling  water  (Model  III).  As  mentioned  above,  a 
special  strut  made  of  pipes  was  used  for  Model  II.  The  dimensions 
and  the  external  geometry  of  the  struts  for  Model  I,  II,  and  III  were 
essentially  the  same. 

Zero  Momentum  Control 

For  the  second  and  third  models,  sufficient  thrust  was  required 
to  make  the  wake  momentumless.  Therefore,  considerable  testirg  was 
needed  to  determine  the  plenum  p -assure  of  the  injected  air  (Model  II) 
or  the  voltage  and  current  of  the  power  supply  driving  the  propeller 
(Model  III)  .corresponding  to  a near  zero  momentum  wake.  Several 
traverses  in  the  left,  right,  and  downward  (opposite  the  strut)  direc- 
tions were  made  in  order  to  find  the  wake  momentum  flux  and  the  sym- 
metry of  the  wake.  The  drag  at  a certain  plenum  pressure  for  Model  II 
or  a certain  voltage  (the  current  was  maintained  constant)  for  Model 
lII  was  calculated  by  the  expression: 

D = 2tt  J*  pu  (u  - Ug)  rdr 


(1) 
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The  pressure  or  the  voltage  that  gave  the  drag  closest  to  zero  was 
chosen  as  the  standard  throughout  five  stations.  This  primary  testing 
was  made  at  station  X/D  ■ 5 because  of  the  low  degree  of  swirl  and  the 
easily  measured  mean  velocity  profiles  at  this  station. 

In  the  second  model,  the  air  injection  was  supplied  by  four  com- 

" • ' , ■' 

pressors  (Ingersoll-Rand  Co.,  400  hp,  ISO  psi,  1600  cfm) . 

A pressure  gauge  was  used  to  control  the  amount  of  the  air  flow.  Care 
has  been  taken  to  keep  the  pressure  constant  because  the  regulator  was 
not  too  efficient.  In  addition  to  the  difficulty  of  controlling 
the  pressure  le*  el,  there  were  also  '»t«er  difficulties  in  contolling 
the  temperature  variation  of  the  air  jet. 

As  it  was  mentioned  previously,  the  propeller  of  the  third  model 
was  driven  by  a motor  which  had  to  be  overloaded  to  render  the  wake 
momentumlcss.  By  using  the  magnetic  pick-up  rpm  display  on  the  oscil- 
loscope, the  speed  of  the  propeller  was  adjusted  by  varying  the  DC 
voltage  supplied  to  the  motor.  Usually  the  propeller  ran  smoothly  and 
did  not  need  any  adjustment;  however,  occasionally  it  would  run 
erratically,  most  probably  due  to  overloading.  In  addition,  it  has 
been  found  that  the  DC  current  was  much  more  suitable  than  AC  current 
for  driving  the  motor,  the  reason  being  that  the  small  fluctuations 
in  rotational  speed  which  were  unavoidable  during  the  AC  operation 
seemed  to  have  a strong  influence  on  the  turbulence  spectrum  produced 
by  the  propeller. 


III.  INSTRUMENTATION 


— 

Comparison  of  Turbulence  Measurement  Techniques 

Basically  there  are  three  broad  classes  of  techniques  for  the 
measurement  of  a fluctuating  field  [20]  - a)  the  u>e  of  a localized 
sending  probe  such  as  a hot  wire  b)  the  use  of  tracers  and  c)  the  use 
of  some  type  of  integrated  measurement,  usually  optical.  The  inte- 
grated methods  cannot  be  used  for  many  applications  because  they 
always  rely  heavily  on  some  additional  assumptions  about  the  statisti- 
cal nature  of  the  field.  The  tracer  techniques  are  very  tedious  and 
time  consuming  because  they  require  many  samples.  For  example,  a 
turbulence  record  of  five  seconds  in  flow  (with  a mean  velocity 
30  m/sec  and  a microscale  of  turbulence  of  0.6  centimeters)  is  equiva- 
lent to  a total  of  25,000  samples.  Therefore,  such  technique  is  not 
convenient  for  a routine  work.  As  for  the  localized  probe  techniques, 
there  are  dark  current  discharge  anemometers,  glow  discharge  anemo- 
meters, streaming  potential  probes  used  as  anemometers,  fluctuating 
total  pressure  probes  and  many  others  besides  hot  wire  anemometers. 

The  hot-wire  anemometer  is  still  considered  as  generally  more  suitable 
than  any  of  the  other  challengers. 

Hot-Wire  Probes 

It  was  mentioned  earlier  that  single-  and  crossed  wires  wore  used 
by  other  researchers.  The  crossed  wire  can  measure  the  turbulence  in 
two  perpendicular  directions  and  the  directional  correlation  (Reynolds 
stress)  at  the  same  time  by  connecting  the  anemometer  with  a correlator 
and  an  RMS  meter  [21].  This  procedure  sounds  handy  and  attractive, 
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but  the  effect  of  asymmetry  is  unavoidable  in  the  crossed-wire  probe 
and  double  sets  of  readings  have  to  be  taken  by  rotating  the  instrument 
180°  about  the  axis  of  flow  direction  between  different  sets.  In 
adlition  to  this,  the  data  reduction  which  depends  on  the  assump- 
tion that  the  turbulent  field  can  be  superposed  linearly  is  not  true 
in  general  [12].  In  this  work,  only  a single  channel  of  the  anemometer 
was  available.  Therefore,  a single-normal  probe  was  used  to  measure 
the  longitudinal  turbulence  intensity,  ana  a single  45°  yawed  probe 
was  used  to  measure  the  Reynolds  stress  in  the  radial  and  tangential 
directions,  and  radial  and  tangential  turbulence  intensity,  by  rotating 
the  wire  180°  about  the  axis  of  flow  direction.  The  probes  used  were 
single  normal  (TSI-1210)  and  single  yawed  probes  (TSI-1213)  shown  in  Figs. 
6 and  7.  The  diameter  and  length  of  the  sensor  were  0.00015  in.  and  0.05 
in.,  respectively.  The  length  of  the  whole  probe  (Figs.  6 and  7)  was 
1.875  in.  The  probes  were  made  of  platinum  plated  tungsten  wires  having 
the  following  characteristics:  1)  a high  temperature  coefficient  of 

resistance  2)  an  electrical  resistance  such  that  it  can  be  easily  heated 
with  an  electrical  current  at  practical  current  and  voltage  levels. 

Anemometer  System 

Turbulence  measurement  data  were  obtained  by  using  a TSI  (Thermo- 
System,  Inc.)  constant  temperature  anerruieter  (Model  1050),  a DC 
Doric  Integrating  Voltmeter  (Model  DS-100) , and  a DISA  type  55D35  RMS 
voltmeter.  The  anemometer  has  very  low  noise  (less  than  0.007% 
equivalent  turbulent  intensity),  high  DC  and  AC  gain  (35,000)  with 
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gain  curve  shaping  control,  high  power  output  (1.5  amps),  and  the 
frequency  response  is  above  500  kHz.  The  output  of  the  anemometer 
was  not  linearized  because  the  mean  velocity  varied  within  a relatively 
small  range.  The  RMS  voltmeter  had  a signals  response  range  of  1 hertz 
to  400  kilohertz,  integrator  time  constants  from  0.1  to  30  seconds,  and 
an  accuracy  of  0.5%  of  the  full  scale  deflection. 

Auxiliary  Equipment 

In  addition  to  the  anemometer  data,  mean  flow  flow  information  [19] 
was  obtained  by  us  ng  a pitot-static  tube  (United  Sensor  USC-A-155)  and 
directional  (United  Sensor  B-llll,  w-187)  and  Kiel  probes. [19]  The  two 
latter  probes  were  used  to  measure  the  inclination  angle  of  the  flow 
and  the  total  pressure,  repectively. 

Automatic  Traverse 

All  probes  were  mounted  on  a traverse  especially  modified  for  this 
investigation  (Fig.  8).  The  probe  could  be  moved  upward  and  downward 
vertically;  left  and  right  horizontally.  The  position  of  each  probe 
was  controlled  by  means  of  variable  speed  electrical  DC  motors  and  was 
indicated  on  a x-y  recorder.  The  recorder  was  calibrated  in  such  a 
way  that  one  inch  movement  of  the  probe  in  either  direction  gave  an 
indication  of  two  inches  on  the  x-y  recording  chart. 


IV  EXPERIMENTAL  PROCEDURE 


Generij. 

In  this  work,  the  main  concern  is  placed  on  the  turbulence  char- 
acteristics - including  axial  (first  subsection),  radial,  and 
tangential  turbulence  intensities,  and  tangential  and  radial  shear 
stress  (second  subsection) . All  the  mean  flow  data  are  given  by 
Swanscn  [19). 

Before  making  the  turbulence  measurements  at  a given  station,  the 
center  of  the  wake  was  determined  by  using  the  pitot-static  tube  and 
making  several  slow  dynamic  pressure  scans  of  the  wake  in  both  the 
horizontal  and  vertical  directions.  The  wake  center  was  identified 
as  a location  of  either  the  maximum  or  minimum  axial  velocity. 


Axial  Turbulence  Intensity 

In  this  experiment,  the  normal  hot  wire  was  used  to  L.aasure  the 
turbulence  in  the  longitudinal  direction.  The  procedure  for  per- 
forming this  part  of  the  measureii.>nts  consisted  of  adjusting  the  hot 
wire  perpendicular  to  the  flow  direction  and  taking  data  at  0.25  inch 
increments  in  the  radial  direction. 


Reynolds  Stress  and  Turbulence  Intensity  in  Radial  and  Tangential 
Directions 

A yawed  hot  wire  was  used  to  measure  the  shear  stresses  and  tur- 
bulence intensity  n the  radial  and  tangential  directions.  The  pro- 
cedure for  performing  this  part  of  the  measurements  consisted  of 
adjusting  the  hot-wire  in  the  vertical  radial  plane  of  the  flow  field 
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(through  the  centerline),  taking  data,  then  turning  the  probe  180° 
about  the  axis  of  the  probe,  and  taking  data  again.  This  was  done  at 
0.25  inch  intervals  across  the  entire  wake.  Next,  the  hot-wire  was 
adjusted  in  the  horiaontal  radial  plane  of  the  flow  field  and  the 
above  procedure  of  taking  data  was  repeated  to  obtain  shear  stress 
and  turbulence  intensity  ir.  the  tangential  direction. 

Calibration  Curves 

The  calibration  of  the  hot-wire  probes  was  made  with  the  aid  of 
a photographic  camera  tripod  mechanism  mounted  in  the  wind  tunnel  test 
section.  The  tripod  was  modified  to  include  a large  protractor  for 
defining  the  inclination  of  the  probe.  The  velocity  sensitivity  was 
determined  by  taking  the  DC  output  at  several  tunnel  q's  and  then 
plotting  the  DC  signal  versus  air  velocity.  The  angular  sensitivity 
of  the  yawed- wire  probes  was  obtained  by  varying  both  the  freestream 
velocity  and  the  angle  of  attack  of  the  probe.  The  angle  of  attack 
was  varied  in  the  range  from  -6°  to  +6*  in  increments  of  2*  A 
sketch  explaining  nomenclature  is  shown  in  Fig.  9.  These  calibrations 
were  done  numerous  times  in  order  to  check  reproducibility  and  to 
obtain  calibration  curves  for  a complete  range  of  temperature. 

Model  Alignment 

Before  the  experiments  were  started,  an  attempt  was  made  to 
align  the  model  with  the  wind-tunnel  stream.  In  the  these  experiments, 
all  three  models  were  aligned  with  the  same  orientation  which  was 
assumed  to  be  parallel  to  the  direction  of  flow  within  an  angle  not 


exceeding  2 minutes.  The  reason  for  the  ^reat  care  needed  to  attain 
an  angle  of  attack  as  close  to  zero  as  possible  is  explained  by 
Buchinskaya  [8] . She  studied  the  structure  of  a vortex  wake  behind 
a thin  bo  iy  of  revolution  at  different  angles  of  attack  and  found  that 
the  structure  of  th;„  vortex  wake  behind  a slender  body  varied  greatly 
with  very  small  changes  in  the  angle  of  attack.  For  our  experiments, 
exact  symmetry  of  the  wake  was  very  difficult  to  achieve.  For  Model  II, 
the  four  brass  interior  vanes  had  some  influence  as  well,  so  they 
were  adjusted  to  have  minimum  effect.  In  addition,  an  imperfect 
geometry  of  the  tear-drop  shape  center  body  could  be  the  most  responsi- 
ble factor  for  a slight  asymmetry  of  the  wake  produced  by  Model  II. 

Temperature  Effects 

Since  the  hot-wire  output  may  be  very  sensitive  to  temperature, 
a particular  care  in  calibration  and  testing  must  be  exercised.  From 
the  calibration  curves  (Fig.  10-12),  it  is  seen  that  1°F  temperature 
difference  in  the  wind  tunnel  caused  an  approximate  0.01  volt  difference 
in  anemometer  DC  output  which  yielded  a 4 ft/sec  difference  in  mean 
velocity  determination.  For  most  of  the  experiments,  the  free  stream 
temperature  variation  remained  within  5°F.  This  5°F  temperature  dif- 
ference implies  as  much  as  a 20  ft/sec  difference  in  the  mean  velocity. 
For  this  reason,  the  mean  velocity  obtained  '..y  the  hot  wire  was  not 
used  directly,  and  the  edge  mean  velocity  obtained  from  the  pitot- 
static.  tube  was  used  to  normalize  the  turbulence  readings.  The  hot 
wires  were  used  only  for  defining  turbulent  properties.  This  is  a 
standard  procedure  for  hot  wire  experiments.  It  is  noticed  that  the 


RMS  output  does  not  vary  much  for  small  temperature  difference  (e.g.^ 
5°F)  . 

Since  Corrsin  and  Uberoi  [24]  found  that  the  initial  jet  tem- 
perature had  a very  great  influence  on  the  temperature  and  velocity 
variation  of  jet  mixing  prollems,  one  must  exercise  great  care  in 
experiments  of  the  type  with  Model  II  here.  Any  temperature  variation 
would  change  the  density  of  the  air  and  hence  the  speed  of  the  jet. 
Therefore,  not  only  the  hot  wire  output  readings  cou’d  be  influenced, 
but  the  speed  of  the  jet  would  also  vary,  so  that  there  was  a possi- 
bility of  having  some  deviation  in  the  mean  velocity  and  turbulence. 

Wire  Contamination 

In  the  tests  on  Model  II,  there  was  some  oil  coming  from  the  com- 
pressors, and  this  might  possibly  contaminate  the  air  stream  and  the 
hot  wire.  Since  the  hot  wire  is  sensitive  to  contamination,  the  data 
obtained  did  not  show  as  good  reproducibility  (0  - 2%  difference  for 
X/D  = 2,  5,  10;  10  - 17%  difference  for  X/D  = 20,  40)  as  for  Models 


V.  DATA  REDUCTION 


The  method  of  obtaining  the  axial  turbulence  intensity  will  be 
described  in  the  first  subsection;  methods  for  obtaining  the  Reynolds 
stresses  and  the  radial,  and  tangential  turbulence  intensities  will 
be  described  in  the  second  subsection.  Basic  to  turbulence  data 
evaluation  were  two  assumptions:  a)  the  flows  approach  those  of 

constant  temperature  and  density  and  b)  the  hot  wires  were  operated 
ideally  by  the  electronic  circuits,  so  that  problem  of  frequency 
response  was  eliminated. 

Axial  Turbulence  Intensity 

The  method  used  for  defining  the  axial  turbulence  intensity  is 
graphical  analysis  [22].  It  is  based  on  the  following  assump-ions : 

1)  the  positive  and  negative  directions  are  the  same  for  velocity  and 
voltage  because  the  hot  wire  anemometer  is  unable  to  distinguish  the 
direction  of  the  flow,  2)  RMS  is  distributed  symmetrically  about  DC. 
At  first,  the  calibration  curve  has  to  be  set  up,  so  that  the  velo- 
city is  known  by  the  corresponding  voltage  output  (See  Fig.  13).  The 
DC  output  gives  the  average  velocity  at  the  "working"  point,  and  the 
RMS  output  gives  the  velocity  fluctuation  around  the  working  point. 
The  longitudinal  turbulence,  >£-2 

= v - v„,  where  v and  v„  are  the 
u B u B 

corresponding  velocities  at  voltage  DC  + 1/2  RMS  and  DC  - 1/2  RMS, 
respectively. 
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Reynolds  Stresses  and  Turhulence  Intensity  in  the  Radi al_ and  Tangential 
Directions 

i)  Radial  Shear  Stress  and  Turbulence  Intensity 

To  calculate  Reynolds  stress  and  turbulence  intensities  in  another 
than  axial  direction,  we  use  the  following  relations: 


i2  = (f)2  u’2  + 2 (|?)(i|£)  (J  ||)2  v'2 


3u'  'u  a*'  “ v ' vu  dip' 


O ry 

e = S Z u'2  + 2S  S u'v ' + S 1 v'2 


U V 


so 


e2  = S 2 u'2  + 2 S .S  . u'v'  + S ? v'2 
+ u+  u+  v+  v+ 


22  = S 2 u'2  + 


e - S u + 2S  S u ' v ' + S 2 v'2 
u-  u-  v-  v- 


(2) 

(3) 

(4) 


where  Sy  = 3E/9U  (velocity  sensitivity),  = 1/v  3E/9^  (angular  sen- 
sitivity) obtained  from  calibration  curves  (see  Experimental  Pro- 
cedures) with  respect  to  velocity  and  angle. 

+ subscripts  represent  the  value  obtained  for  the  slanted  probe 
while  the  long  prong  is  on  the  top  (wire  inclined  +45°) 

- subscripts  represent  the  value  obtained  for  the  slanted  probe 
while  the  short  prong  is  on  the  top  (wire  inclined  -45°) 

The  two  unknowns  (\Zv"'2,  u'v')  can  be  obtained  from  the  equa- 
tions 3 and  4. 

1 i ) The  Tangential  Shear  Stress  and  the  Turbulence  Intensity 

I hose  two  quantities  are  obtained  by  the  same  method  as  above 
except  that  the  + and  - subscripts  represent  the  values  obtained  with 


respectively,  and  v in  the  equations  is  replaced  fcy  w. 


I 1 1 


VI  RESULTS  AND  DISCUSSION 


As  described  in  the  foregoing,  three  different,  models  were  investi- 
gated. The  first  model  is  a streamlined  body,  and  the  flow  field 
behind  this  body  is  simpler  than  the  other  two.  The  flow  has  lower 
turbulence  intensity  and  Reynolds  stress,  and  the  wake  is  produced  by 
the  body  only.  The  wake  produced  by  the  second  model  is  a mining  of 
the  air  jet  and  the  wake.  This  involves  the  characteristics  of  the  air 
jet  and  the  geometry  itself.  The  wake  produced  by  the  third  model  is 
more  complicated  than  those  of  the  other  two,  because  the  propeller 
produces  velocities  in  both  the  tangential  and  radial  directions,  and 
the  problem  involves  the  turbulence  characteristics  of  both  the  pro- 
peller and  the  body  wake. 

Since  all  these  three  models  are  essentially  axisymmetric , the 
wakes  induced  should  be  al  so  axisymmetric  provided  that  the  alignment 
and  the  construction  are  perfect.  Several  mean  flow  measurements  were 
made  to  test  the  symmetry  of  the  wakes  and  it  was  found  that,  in  general, 
the  mean  flow  symmetry  was  satisfactory  [19].  In  addition,  turbulence 
measurements  in  the  axial  direction  at  various  stations  were  made 
to  assess  the  symmetry  of  the  turbulence  distribution  in  the  wake  of 
Models  I and  II  (see  Appendix  A).  In  the  case  of  Model  I,  the  "tur- 
bulence symmetry"  was  judged  satisfactory  (cf.  Table  A-l) . For  the 
case  of  Model  II,  the  wake  was  fairly  close  to  symmetry  except  for  the 
outer  portion  of  the  wake  (r  = 2-3  in.)  at  near  stations  (X/l)  = 2 and 
rj).  Asymmetry  at  the  locations  just  mentioned  might  be  caused  by 
slight  misalignment  and  asymmetry  of  the  four  vanes  in  the  peripheral 
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slot  of  Model  II.  The  downward  direction  (opposite  to  the  strut)  was 
selected  for  detailed  measurements  and  w?  . taken  as  representative  of 
other  directions. 

In  this  work,  the  turbulence  intensities  and  Reynolds  stresses 
are  the  main  subject  of  interest  and  will  be  discussed  separately  for 
the  different  models.  All  the  mean  flow  data  are  presented  and  dis- 
cussed in  Ref.  19. 

Model  I:  Streamlined  Drag  Body 

Model  I is  the  only  one  of  the  three  models,  which  does  not  pro- 
duce a momentumless  wake,  and  as  such  it  serves  as  a basis  for  com- 
parison of  the  self-propelled  cases.  The  turbulence  intensity  profiles 
of  the  wake  in  the  three  perpendicular  directions  are  plotted  in 

Figs.  14,  15,  and  16.  The  Reynolds  stress  profiles  are  plotted  in 

* 

Figs.  17  and  18.  From  these  plots,  it  is  apparent  that  the  wake  spreads 
as  it  moves  downstream.  The  axial  turbulence  intensity  profiles  (Fig. 

14)  show  that  the  peak  value  increases  over  a distance  ranging  from 
12  inches  (X/l)  = 2)  behind  the  body  tail  to  60  inches  (X/D  = 10)  and 
then  decreases  as  it  moves  further  downstream.  The  other  two  direc- 
tional turbulence  intensities  exhibit  similar  trends;  the  peak  tur- 
bulence intensity  increases  along  a distance  extending  from  12  inches 
behind  the  body  tail  to  30  inches  (X/D  = 6)  and  then  drops  as  the 
wake  moves  further  downstream.  The  major  reason  for  such  a behavior 
lies  in  the  convergence  of  the  wake  near  the  body  tail.  By  comparing 
the  three  figures  (14,  15  and  16),  it  can  be  seen  that  the  longitudinal 
turbulence  intensities  are  higher  than  those  in  the  other  two  directions. 

* Alt  data  plotted  are  corrected  for  background  turbulence. 
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The  radial  Reynolds  stress  is  twice  as  high  as  that  in  the  tan- 
gential direction  (Fig.  18).  The  presence  of  the  tangential  shear 

stress  may  indicate  some  asymmetry  of  the  wake  due  to  an  imperfect 

2 

alignment.  The  peak  radial  Reynolds  stress  (u'v'/u  ) increases 

e max 

up  to  station  X/D  = 5 and  then  decreases  downstream. 

Now,  the  results  will  be  compared  with  those  obtained  by  Chevray 
and  Bukreev  [10,  11].  When  comparing,  one  has  to  realize  that  the  body 
used  by  Chevray  is  somewhat  similar  to  this  model  (Model  I) , the  main 
diffe'ence  being  that  the  tail  of  this  model  is  sharp  while  Chevray's 
is  blunt  and  the  ratio  of  the  length  to  the  maximum  diameter  of  his 
model  is  6:1,  while  for  the  Model  I this  ratio  is  12:1.  The  body 
studied  by  Bukreev  is  sharp  tailed.  Table  II  and  Fig.  19  present  the 
representative  values  obtained  in  the  three  experiments.  It  can  be 
seen  that  the  results  being  compared  are  close  to  each  other  in 
magnitude,  and  that  the  longitudinal  turbulence  intensities  arc  higher 
than  those  in  the  other  two  directions.  In  the  results  obtained  by 
Chevray,  the  peak  turbulence  intensities  in  the  radial  direction  and 
the  peak  Reynolds  stresses  decrease  starting  from  X/D  = 9.  Such 
phenomena  occur  due  to  a complex  vortex  system  formed  in  the 
boundary  layer  of  the  body,  and  converging  with  the  flow  direction 
in  the  near  wake  [13].  The  results  obtained  in  this  work  are  some- 
what different  which  can  be  explained  by  the  fact  that  the  tail  of  the 
two  models  are  different  in  shape,  so  that  the  vortex  systems  converge 
< 1 i f ferently . 

The  similarity  profile  for  the  axial  turbulence  intensity  is 


4 


i 

1 


, 


i 


shown  In  Fig.  20.  This  profile  is  formed  by  normalizing  the 


Representative  Turbulent  Quantities  for  the  Present  Drag  Body  (Model  I)  and  for  the  Models  Studied 
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turbulence  intensity  and  the  radial  distance  with  the  values  of 


(//u'2)  and  R*,  respectively,  where  (Vu'2)  is  the  maximum  axial 
max  max  _ 

turbulence  and  R*  is  the  r-distance  corresponding  to J u'2=  ~(//u'2) 

2.  max 

The  different  distributions  form  a unique  curve  in  the  outside  portion 
of  the  wake;  the  scattered  data  at  the  center  of  the  wake  confirm 
that  the  turbulence  fluctuations  are  not  self-preserved.  The  simi- 
larity profiles  for  radial  shear  stress  are  shown  in  Fig.  21.  The 


similarity  profile  is  formed  by  the  scaling  factor  (u'v')  and  R*. 

max 

The  different  distributions  of  these  profiles  form  a unique  curve  in 
both  the  inner  and  the  outer  regions  of  L.*e  wake  with  some  scatter. 
All  these  results  are  similar  to  the  results  of  Chevray. 


Model  II  : Jet-Propelled  Bc^y 

The  turbulence  intensity  profiles  are  plotted  in  Figs.  22,  23, 
and  24.  These  figures  show  that  the  profiles  are  very  different  from 
those  for  the  streamlined  body  both  in  their  distributions  and  t’  air 
magnitudes  for  all  representative  quantities.  The  longitudinal  tur- 
bulence intensity  profiles  have  some  peculiar  phenomenora  near  the 
center  of  the  wake  and  there  are  dips  near  r/R  = 0.25  at  the  stations 
close  to  the  tail  (X/l)  = 2 and  X/L  = 5).  The  radial  and  tangential 
turbulence  intensity  profiles  are  not  smooth  in  the  inner  portion  of 
the  wake.  The  reason  for  this  may  be  related  to  the  character  of  the 
air  jet.  The  jet  passed  screens,  converged  along  the  tear-drop  shaped 
renter  body  an  I then  mixed  toget  her  - which  is  a complex  system, 

I mated  in  tin  Inner  portion  of  the  wake . The  outer  portion  of  the 
wake  is  produced  primarily  by  the  main  body  and  is  smoother.  The 
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Reynolds  stress  profiles  (Figs.  25  and  26)  show  that  there  is  a very 
hig  f ip  after  the  first  peak  but  no  such  dip  appears  in  the  stream- 
lined body  data.  Note  that  the  sign  of  the  shear  stress  is  opposite 
to  that  of  the  streamlined  body.  It  is  obvious  that  installing  the 
air  jet  can  change  the  whole  flow  field. 

For  Model  II, the  wake  is  wider  than  that  of  the  streamlined  drag 
body.  The  peak  turbulence  intensities  in  the  longitudinal,  radial 
and  tangential  directions,  and  the  peak  Reynolds  stress  decrease, 
as  the  wake  moves  downstream.  The  turbulence  intensity  in  the  longi- 
tudinal direction  is  higher  than  that  in  either  of  the  other  two 
directions  (see  Figs.  22-24). 

Since  there  are  no  published  experimental  data  on  a slender, 
streamlined  body  with  momentumless  wake  available,  the  work  used  for 
comparison  here  is  a study  of  a circular  disk  with  air  injection,  which 
is  neither  slender  no  streamlined  (Naudascher,  Ref.  12).  Table  III  and 
Fig.  27  show  that  the  peak  turbulence  intensity  and  Reynolds  stress  for 
Model  II  are  lower  than  for  Naudascher's  models.  The  longitudinal  tur- 
bulence intensities  are  higher  than  the  radial  and  tangential  ones,  but 
the  peak  radial  and  tangential  turbulence  intensities  measured  here 
are  close  to  the  peak  axial  turbulence  intensities  for  Naudascher's  model. 
The  decay  rates  of  the  turbulence  intensity  of  Model  II  are  much  slower 
than  in  Naudascher's  experiments.  The  distribution  of  longitudinal  tur- 
bulence intensities  and  the  shape  of  Reynolds  stress  profiles  are  suf- 
ficiently similar  to  each  other.  It  is  noticed  that  the  dip  is  higher 
than  the  peak  in  the  Reynolds  stress  profiles  for  this  model  but  the  re- 
vet-,e  is  true  for  Naudascher's  model.  The  differences  may  be  due  to 
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the  basic  differences  in  the  model  geometries  and,  hence,  wakes 
produced  by  Model  II  and  Naudascher's  model. 

The  self-similarity  profiles  for  axial  turbulence  intensities  and 
radial  shear  stress  are  shown  in  Figs.  28  and  29.  The  axial  turbulence 
intensity  distributions  can  form  a definite  curve  in  the  outside  portion 
of  the  wake  although  there  is  some  scatter  at  the  station  X/D  = 40. 

The  reason  for  the  scatter  at  this  station  is  that  its  turbulence 
intensity  is  very  low  as  compared  with  that  at  other  stations,  and  the 

H—  9 

experimental  error  is  amplified  by  the  factor  (v  u'  ) . The  scatter 

max 

in  the  inner  portion  of  the  wake  is  caused  by  the  characteristics  of 
the  air  jet.  It  is  seen  that  the  radial  shear  stress  distribution  does 
not  form  a good  definite  curve  in  the  outside  portion  of  the  wake.  It 
seems  as  if  they  form  two  curves  - one  for  stations  X/D  = 2 and  5,  and 
the  other  for  X/D  = 10,  20,  and  40.  The  reason  for  the  observed  ir- 
regularities may  be  that  the  air  jet  has  a larger  influence  over  near 
stations  than  the  free  stream;  at  the  farther  stations,  the  air  jet  and 

the  free  stream  flow  mix  more  and  both  are  important  to  the  outside  portion 
of  the  wake. 

Model  III:  Propeller-Driven  Body 

The  turbulence  intensity  profiles  in  the  longitudinal,  radial  and 
tangential  directions  are  plotted  in  Figs.  30-32.  The  radius  of  the 
wake  is  much  wider  than  that  for  the  other  two  models.  The  longitu- 
dinal turbulence  intensity  profiles  have  more  than  one  peak  at 
stat  Ions  X/D  = 2,  5,  and  10.  The  radial  and  tangential  turbulence 


Representative  Turbulent  Quantities  for  the  Present  Propeller-Driven  Body  (Model  III)  and  for 

the  Mode]  Studied  by  Gran  [14] 
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intenFity  profiles  each  have  a second  peak  after  the  first  one  and 
the  second  peaks  are  lower  at  downstream  stations.  Moreover,  the 
height  of  the  first  peak  decreases  as  the  wake  moves  downstream. 

Another  significant  point  to  be  noticed  is  that  the  turbulence  inten- 
sity profiles  are  quite  flat  and  widely  spread  after  X/D  «-  10.  The 
Reynolds  stress  profiles  (Fig.  33)  in  the  radial  direction  each  have 
a peak  and  then  dip,  opposite  to  that  observed  for  the  second  model 

I 

(see  Figs.  25  and  33).  The  magnitude  of  the  peak  is  bigger  than  that 
of  the  dip  but  all  of  them  decay  very  rapidly  downstream  after  the 
station  X/D  = 5.  Domstream  of  X/D  = 10,  there  is  no  dip,  the  profiles 
become  very  flat  and  the  values  are  very  low  as  compared  with  those  at 
stations  near  the  boby.  Reynolds  stress  values  in  the  tangential 
direction  (Fig.  34)  are  bigger  than  for  the  other  two  models.  We 


observed  an  even  higher  tangential  shear  stress  than  that  in  the  radial 
direction  at  the  stations  near  the  tail  of  the  body.  This  is  due  to 
the  presence  of  more  swirling  and  non-negligible  tangential  velocity 
in  the  flowfield  due  to  the  propeller.  In  addition,  the  decay  rates 
of  peak  turbulence  intensities  and  Reynolds  stress  for  this  model  are 
much  higher  than  the  corresponding  values  for  the  other  two  models 
(Figs.  36-40),  and  the  longitudinal  turbulence  intensities  (Figs.  36- 
38,  and  31-33)  are  lower  than  those  in  the  other  two  directions  for 
this  model,  while  for  the  other  two  models  the  longitudinal  turbulence 
intensity  is  higher  than  that  in  the  other  two  directions.  It  is 
noted  that  there  is  some  peculiar  behavior  (Fig.  34)  near  the  center 
of  the  wake,  possibly  due  to  the  effect  of  the  tail  of  the  body  and 


the  hub  of  the  propeller.  It  can  be  concluded  that  all  the  turbulence 
characteristics  of  the  wake  of  Model  III  are  very  different  from  the 
wake  turbulence  characteristics  of  Models  I and  II  both  in  magnitude 
and  profile  distributions  of  turbulence  intensity  and  shear  stress. 

The  turbulence  data  for  Model  III  will  be  compared  with  the 

results  obtained  by  Gran  [14].  Only  ( ) and  OTVVU  2) 

e max  e 'max 

are  listed  in  Table  IV  and  plotted  in  Fig.  35,  because  Gran  measured 
only  axial  turbulence  intensity  and  ;adial  shear  stress.  This  com- 
parison shows  that,  in  terms  of  the  order  of  magnitude  of  the  turbulence 
intensity  and  the  radial  shear  stress,  our  experimental  results  ar; 
generally  in  rough  agreement  with  those  of  Gran.  Here  the  axial  Lur- 
bulence  intensity  of  Model  III  exhibits  more  than  one  peak  at  the 
stations  near  the  body  (X/D  = 2 and  5),  and  the  corresponding  profiles 
are  very  flat  at  the  stations  further  downstream  (X/D  = 10,  20,  and 
40).  The  axial  turbulence  intensity  profiles  obtained  by  Gran  did  noc 
have  a second  peak  but  they  did  have  some  small  ripples.  It  is  seen 
from  Fig.  35  that  the  peak  radial  shear  stresses  for  the  two  models 
are  close  in  order  of  magnitude  but  the  axial  turbulence  intensity  for 
Gran's  model  is  higher  than  that  for  Model  III  although  the  decay 
rates  are  similar.  The  differences  shown  might  be  caused  by  the  dif- 
ferent characteristics  and  performance  of  the  propellers. 

Comparison  of  the  Results  for  Models  I,  II,  and  III 

Wien  the  maximum  values  of  axial  turbulence  intensities  are  com- 
paied  for  the  three  models  (Fig.  36),  it  can  be  seen  that  the  air 
in  lection  model  has  the  highest  value.  The  propeller-driven  model 
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produces  higher  axial  turbulence  intensity  than  the  drag  body  at 
the  nearest  station,  but  the  axial  turbulence  intensity  for  Model  III 
decreases  very  fast  after  station  X/D  = 10  while  this  quantity  for 
Model  I increases  a little  up  to  station  X/D  = 10  and  then  decreases 
at  a similar  rate  as  for  Model  II.  Therefore,  the  propeller-driven 
body  generates  the  lowest  axial  turbulence  intensity  after  station 
X/D  = 10. 

From  Figs.  37  and  38,  it  is  seen  that  the  propeller-driven  body 
produces  the  highest  radial  and  tangential  turbulence  intensities 
and  their  rate  of  decrease  is  faster  than  for  the  drag  body.  Con- 
sequently, the  turbulence  intensities  in  radial  and  tangential  direc- 
tions are  close  for  these  two  models  at  stations  X/D  » 20  and  40.  The 
radial  and  tangential  turbulence  intensities  for  Model  II  at  near 
stations  (X/D  = 2 and  5)  have  values  in  between  those  for  two  other 
two  models  but  they  decrease  slowly. 

From  Fig.  39,  it  is  seen  that  for  the  self-propelled  bodies  the 
initial  magnitude  of  the  radial  shear  stress  is  higher  than  for  the 
drag  body.  For  the  propeller-driven  body  the  radial  shear  stress 
decreases  very  fast,  so  that,  at  downstream  locations,  it  becomes 
even  lower  than  that  for  the  drag  body.  The  radial  shear  stress  for 
the  air-injection  body  decreases  somewhat  faster  than  that  for  Model  I, 
but  is  is  still  higher  than  that  for  Model  I up  to  the  station  X/D  = 
40. 

From  Fig.  40,  it  can  be  seen  that  tangential  shear  stresses  in 
the  wakes  of  Models  I and  II  are  low  but  not  negligible.  This  may 
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be  related  to  some  degree  of  non-uniformity  and  imperfect  symmetry 
in  the  flow  field  behind  these  models.  The  decay  rates  of  the  tan- 
gential shear  stresses  for  the  three  models  are  similar  at  stations 
downstream  of  X/D  = 10. 

In  Fig.  41,  R is  the  distance  from  the  center  of  the  wake  to 
the  location  where  the  axial  turbulence  approaches  free  stream  tur- 
bulence. This  figure  shows  that  the  propeller-driven  body  produces 
the  widest  wake  throughout  the  five  measurement  stations,  and  the 
width  of  the  wake  of  the  streamlined  drag  body  is  the  smallest  one 
except  for  the  locations  X/D  > 20. 

Tho  Figs.  36  through  40  suggest  that  air  injection  increases  the 
turbulence  intensities  in  all  directions;  it  also  increases  the  radial 
and  tangential  shear  stresses.  The  propeller  decreases  the  axial  tur- 
bulence intensity,  but  increases  the  turbuience  intensity  in  the 
radial  and  tangential  directions  and  the  '/angential  shear  stress. 

Momentum  Analysis 

In  order  to  document  the  wake  conditions  and  obtain  an  insight,  into 
the  flow  structure  and  the  mechanism  of  momentum  and  energy  transformation, 
momentum  and  energy  balance  analyses  of  the  data  are  presented.  The 
momentum  analysis  can  be  used  as  a measure  to  determine  whether  the  ex- 
perimental wake  is  the  desired  wake  for  the  models  tested.  A useful 
form  of  the  momentum  equation  is  derived  by  integrating  the  first 
equations  of  Reynolds  over  a cylindrical  control  surface  and 


(1 


(5) 


The  first  term  on  the  right  side  of  Eq . (5)  represents  the  difference 
in  the  mean  momentum  flux  between  the  final  and  initial  sections;  the 
second  term  represents  the  net  flux  of  momentum  through  the  periphery 
of  the  region  under  consideration  and  the  third  term  represents  the 
effect  of  the  pressure  difference  between  two  sections;  the  last  term 
represents  the  difference  in  the  turbulence  momentum  flux. 

We  define  - u;  then  Eq . (5)  becomes 
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Hence 

-cD  - 4 <I2  - Ti  - V + 2Ip 

The  values  of  integrals  are  tabulated  in  Tables  V-VII  and  plotted 
in  Figs.  42-44. 

i)  Model  I:  ST'amlined  Drag  Body 

The  results  are  shown  in  Table  V and  Fig.  42.  Theoretically,  the 
drag  coefficient  at  various  stations  should  be  exactly  equal,  but  in 
these  experiments  there  are  deviations  up  to  15%  from  the  final  drag, 
coefficient,  = 0.104.  However,  it  can  be  seen  that  the  drag 
approaches  a constant  value  for  the  downstream  stations.  The  non- 
negligible  tangential  force  may  be,  at  least  in  part,  responsible  for 
discrepancies  at  the  near  stations. 

ii)  Model  II:  Jet-Propelled  Body 

The  results  are  shown  in  Table  VI  and  Fig.  43.  For  a zero  momen- 
tum wake,  the  drag  should  be  zero.  At  stations  other  than  X/D  = 2,  the 
drag  coefficient  is  very  close  to  zero  while  at  station  X/D  = 2 it  has 
a value  of  0.005.  Such  discrepancy  may  be  related  to  a complicated 
flow-field  near  the  tail  of  the  body,  and  difficulty  of  measurements 
at  this  station. 

iii)  Model  III:  Propeller-Driven  Body 

The  results  are  shown  in  Table  VII  and  Fig.  44.  From  Table  VII, 
it  is  seen  that  the  drag  coetficient  is  close  to  zero  except  for  the 
station  X/D  = 2,  where  = 0.021.  Such  discrepancy  at  X/D  = 2 may 
be  caused  by  the  large  tangential  velocities  in  the  flow  field  near  the 
t iii  of  the  body  and  again  difficulties  of  measurements  at  this  station. 

1 rom  Tallies  V-VII,  it  can  be  seen  that  the  turbulence  momentum  and 
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TABLE  V 

Momentum  Integrals  : Model  I 


X/D 

X1 

X2 

IT 

1 = 
P D 

4drI2-IT)-2Ip 

2 

0.0299 

0.00725 

0 00034 

0.000152 

0.0869 

5 

0.0252 

0.00508 

0.000352 

0.000547 

0.07798 

10 

0.0275 

0.00416 

0.000576 

0.0023 

0.08644 

20 

0.0309 

0.00263 

0.000875 

0.0027 

0.1042 

40 

0.0288 

0.00111 

0.00908 

0.0015 

0.1040 

TABLE  VI 

Momentum  Integrals:  Model  II 

X/D 

X1 

X2 

XT 

I 

P 

cd  ' 4 (IrVV-2Ip 

2 

0.00585 

0.00186 

0.0024 

0.000712 

0.0049 

5 

0.00333 

0.001327 

0.00149 

0.000896 

0.00026 

10 

0.00254 

0.00302 

0.00102 

0.002712 

0.00108 

20 

0.00265 

0.000115 

0.00129 

0.002240 

0.00137 

40 

0.00017 

0.000002 

0.0004615 

0.000235 

0.001644 

TABLE  VII 

Momencura  Integrals:  Model  III 

X/D 

X1 

X2 

XT 

L CD- 

4 “rW-^p 

2 

0.00387 

0.00425 

0.001174 

0.00748 

-0.02117 

5 

0.00511 

0.00152 

0.0089 

0.00636 

-0.00192 

10 

0.00478 

0.000408 

0.000885 

0.00807 

-0.002192 

20 

0.00470 

0.00007 

0.00032 

0.00763 

0.00254 

40 

0.00382 

0.000003 

0-00009 

0.00585 

0.0032 
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pressure  contribution  cannot  be  neglected,  and  chat  the  air-jet  speed 
in  the  second  model  and  the  rpm  of  the  propeller  in  the  third  model  are 
sufficient  to  create  nearly  momentumless  wake.  From  Figs.  42-44,  it  is 
seen  that  there  are  changes  in  all  values  of  momentum  integrals  near 
the  tail  of  the  three  bodies.  This  indicates  the  direct  contribution 
of  the  body  itself.  The  values  of  the  momentum  integrals  approach  con- 
stant values  far  downstream,  except  for  some  deviations  for  the  case 
of  Model  II.  The  drag  coefficients  for  Model  II  and  Model  III  are 
fairly  close  to  zero.  The  momentum  change  due  to  the  mean  flow  plays  an 
important  role  in  the  momentum  analysis  (see  Fig.  42-44),  so  that  the 
accuracy  of  mean  flow  measurements  can  be  assessed  by  the  momentum 
inn  lysis. 

Energy  Balance  Analysis 

Assuming  that  the  Reynolds  number  is  sufficiently  high  so  that  the 
work  done  by  viscous  stresses  is  negligible,  one  can  write  the  energy 
equation  for  a steady,  axisymmetric  mean  flow  of  an  incompressible 
fluid  without  a spiral  component  [26]: 

2/fc,l^i)-1]td5  + 2/  8-  r 

, C , u u'2  v u'v'.r  r C u r r 

+ V (v  + r rr) « d K + 1/ (1  - K d i 

e uo  e ue  ^ e 

r r\~ ~»  *u,u  .v,u  , y ,u  ,u 

+ 4 / / ii-v  (-LL* K + 9 1 e)  + ±J.  e 

JJ  Lue2  Wr  9x/R  ; 3x/r 

, v:.2  <^V/1  e w' 2 V^Ue|  r_  . r_  , x _ _ . 

U 2 9r/R  iT1  r/R  J R d R d R ° (7) 
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The  first  and  fourth  terms  come  from  the  energy  equation  of  mean 
motion.  The  second,  third  and  fifth  terms  represent  the  energy  con- 
tributions due  to  pressure  differences,  turbulence  and  shear  stresses. 
For  convenience,  the  following  symbols  are  used: 

2.  2 


‘■7 

/ 


,2 

c i y__  y__  l a l 

2 / u „ 2 R d R 

e LI 


/ 


c = i L L j I 

S3  _/  U R d R 
e 


f f u'v'  9U/Ue  r r x 

4 JJ  Ug2  3r/R  R R R 

f f L A L A * - fc  2 

5 JJ  Ug2  3x/R  R d R d R J S6  d P 

/ v_  u'v'  r , £ 

1 J Ue  Ue2  R R 

= ffn'v'  3V/Ue  r r x 

8 yy  ue2  8x/R  R R R 


10 


Ue2 

3x/R 

’ v'  2 

.V/U 
3 / e 

».2 

3r/R 

/U 

v/  e 

Ue2 

r/R 
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Eqn.  (7)  can  be  expressed  as 


4S10  + 2S 


11  + 2S1  + 4S2  + 4S7  “ 2S3  + AS  + AS.  + 4S,  + 


AS, 


(8) 


The  values  of  energy  integrals  are  tabulated  In  Table  X-XII  and 

plotted  in  Figs.  45-47  (A  is  the  deviation  of  the  sun,  of  integrals 
in  Eqn . (7)  from  zero). 

i)  Model  I:  Streamlined  Drag  Body 

For  the  streamlined  body,  the  radial  velocity  v Is  assumed  to  be 
negligible  as  compared  with  the  axial  velocity  at  the  stations  X/D  - 2 
through  40;  therefore,  only  six  terms  are  considered  (Sj , s;,  Sj,  S(, 


s5,  and  su).  Table  VIII  shows  that  there  Is  only  a small  deviation 


(A)  of  the  experimental  result  as  compared  with  the  energy  balance 
equation.  Moreover,  It  shows  that  the  turbulence  contribution  is  very 
small  as  compared  to  the  mean  velocity  contribution,  and  the  energy 
balance  in  this  case  is  not  sensitive  enough  to  detect  errors  In 
turbulence  measurements. 

ii)  Model  T I : Jet-Propelled  Body 

Since  the  flow  field  of  the  air-injection  model  is  more  com- 
pJirated  than  that  of  the  streamlined  body,  the  radial  velocity  may 

'K,V'  t0  bC  laken  int°  acCo,,nt*  However,  the  radial  velocity  was  not 
tied  Irom  the  measurements  and,  therefore,  some  estimates  are 
needed.  By  the  continuity  equation  in  the  cylindrical  coordinate 


I 


d 


TABLE  VIII 

Energy  Integrals:  Model  I 


X/D 


2S 


11 


2 

5 

10 

20 

40 


0.00315 

-0.00225 

0.0009 

0.00145 

0.0035 


0.47583 

0.86841 

0.86573 

1.0973 

1.3622 


0.44215 

0.83812 

0.82882 

1.0485 

1.3365 


-0.796x10 


-5 


-0.188x10 


-4 


-0.319x10 


-4 


0.00030  -0.000192  -0.06042 

0.06576 


-0.719x10 


-4 


0.00038  -0.00055 
0.00051  -0.00115 


-0.135x10 


-3 


0.00083  -0.00146 
0.00066  -0.00211 


-0.07458 

-0.09867 

-0.0502 
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TABLE  IX 


Radial 

Velocity  of 

Estimates  for 

Model  II 

v/uSjC/D 

? 

5 

10 

20 

40 

r/R>v\ 

0.1 

0.00115 

0.000165 

o.or 

-0.000056 

-0.000032 

0.2 

0.000542 

-0.000110 

-0.000165 

-0.000052 

-0.000062 

0.3 

0.000251 

-0.000108 

-0.00030 

-0.000078 

-0.000085 

C.4 

-0.000347 

-0.000615 

-0.000325 

-0.000087 

-0.000082 

0.5 

-0.000642 

-0.00070 

-0.000355 

-0.000091 

-0.000076 

0.6 

-0.000415 

-0.000515 

-0.000265 

-0.000058 

-0.000067 

0.7 

0.000201 

-0.000214 

-0.000097 

-0.000002 

-0.000061 

0.8 

0.000738 

-0.000055 

0.00001 

0.0000475 

-0.000051 

0.9 

0.00080 

0.0001825 

0.0001411 

0.0000918 

-0.000043 

1.0 

0.0007817 

0.0002117 

0.0001269 

0.0000826 

0.0000387 

TABLE  X 

Energy  Integrals:  Model  II 


X/D 

S11 

S1 

S2 

S4 

S3 

A 

2 

0.345x10'^ 

0.78963 

0.00218 

-0.000932 

0.78353 

0.0192 

5 

0.15  xl0-‘ 

1 4030 

0.0017298 

-0.00195 

1.3922 

0.0204 

10 

0.25  xlO  ^ 

1.3974 

0.001362 

-0.00263 

1.3914 

0.0064 

20 

0.678x10 

2.0028 

0.001854 

-0.00288 

2.00844 

0.00721 

40 

0.410x10 

2.0005 

0.000569 

-0.00308 

^ . 0002 

0.0086 

40 


system  written  in  the  nondimens ional  form: 

3U/Ue  2 3_  , v_  r, 

3x/D  “ ' r/R  3r/R  VUg  R 

v_  11  f1  3U/9Ue  r r 

U " r/R  2 J 3x/D  R R 

e o 


(9) 


the  radial  velocity  v can  be  determined.  This  is  done  by  plotting  the 
axial  mean  velocity  vs  x/D  at  a fixed  r/R  and  calculating  and  then 
integrating  it.  The  results  are  listed  in  Table  IX.  This  table  shows 
that  the  radial  mean  velocity  is  very  small  as  compared  with  the  axial 
velocity.  Therefore,  all  the  terms  involving  v/Ue  in  the  energy 
balance  equation  can  be  neglected  so  that  only  six  terms  (S^,  S2,  S^, 

S , S , S,,)  are  considered.  Of  these,  S is  roughly  two  orders  of 
magnitude  smaller  than  the  remaining  terms  and,  therefore,  it  can  be 
neglected.  The  value  of  energy  integrals  for  Model  II  are  listed  in 
Table  X and  plotted  in  Fig.  46.  From  Table  X,  it  is  seen  that  the 
deviation  of  the  energy  Dalance  equation  is  even  smaller  than  the 
deviation  shown  for  Model  I.  Also,  the  turbulence  contributions  are 
larger  while  the  pressure  contributions  are  smaller  than  those  for 
Model  I. 

iii)  Model  III:  Propeller-Driven  Body 

The  flow  field  of  this  model  is  not  simple,  hence  all  the  terms 
of  equation  (7)  should  be  evaluated  as  for  Model  II.  The  radial  mean 
velocity  was  esimated  by  the  same  method  as  for  Model  II,  since  the 
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degree  of  swirling  is  not  high  (the  yaw  angle  is  less  than  seven 
degrees  at  the  station  X/D  = 2)119].  The  radial  velocity  estimated  as 
above  is  listed  in  Table  XI.  From  this  table,  it  Is  seen  that  the 
radial  velocity  is  larger  than  that  in  the  air-inject ton  case.  However, 
the  values  are  still  very  low  as  compared  with  the  axial  velocity, 
therefore,  it  is  still  safe  to  neglect  all  the  terms  including  v/u^  in 
the  energy  balance  equation.  The  term 


u ,U 
/ e 

r/R 


is  negligible  also.  The  values  of  energy  integrals  for  Model  III  are 
listed  in  Table  XII  and  plotted  in  Fig.  47.  Table  XII  shows  that  the 
deviations  of  experimental  results  are  small  as  compared  to  those  for 
Model  I but  they  are  slightly  larger  than  or  equal  to  the  deviations 
estimated  for  Model  II.  The  pressure  and  turbulence  contributions  are 
larger  than  those  for  Model  I.  The  deviation  A is  the  largest  at 
station  X/D  = 2 where  swirl  is  the  largest . 


I. 
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TABLE  XI 

Radial  Velocity  Estimates  for  Model  III 


u.l 

0.001415 

0.00095 

0.2 

0.002245 

0.0016675 

0.3 

0.00212 

0.0018575 

0.4 

0.0034? 

0. 0017425 

0.5 

0.002822 

0.001443 

0.6 

0.002822 

0.0011566 

0.7 

0.0021128 

0.000666 

0.8 

0.0015299 

0.0002455 

0.9 

0.0010 

TABLE  XII 


Energy 

Integrals : 

Model  III 

X/D 

S11 

s 

S1 

S2 

S4 

A 

2 

-0.0052 

2.2/261 

2.7460 

0.00127 

-0.00088 

-0.03096 

5 

-0.0036 

2.0051 

2.0198 

0.00189 

-0.00195 

0.022 

10 

0.0016 

2.0048 

2.0156 

0.00169 

-0.00277 

-0.02048 

20 

C. 00351 

3.5508 

3.5416 

0.00045 

-0.00307 

-0.0218 

40 

0.0039 

3.5517 

3.5442 

0.00027 

-0.00308 

-0.01844 

VII.  CONCLUSIONS 


In  this  work,  the  turbulence  properties  were  measured  in  the  wake 
behind  three  separate,  but  related  bodies:  Streamlined  Drag  Body, 

Jet-Propelled  Body  and  Propeller-Driven  Body;  all  with  the  same  fore- 
body shape.  The  wake  produced  by  the  latter  two  models  was  momentum- 
less.  The  turbulence  properties  included  the  axial,  radial,  and 
tangential  turbulence  intensities  and  the  radial  and  tangential  shear 
stresses,  and  the  measurements  were  made  at  five  stations:  X/D  = 2, 

5,  10,  20  and  40.  After  comparing  the  present  results  of  the  three 
models  with  each  other  and  with  the  experimental  results  by  other 
researchers,  general  conclusions  may  be  drawn. 

First,  since  the  flow  fields  behind  the  tested  bodies  are  very 
different,  the  turbulent  characteristics  of  the  wake  of  the  three 
models  are  very  different  from  one  another.  The  size  of  the  wake  and 
the  rates  at  which  they  decay  are  different.  The  wake  of  Model  I is 
narrower  than  that  of  Model  II  which  in  turn  is  narrower  than  that  of 
Model  III.  The  axial  turbulence  intensity  profiles  measured  in  the 
radial  direction  are  smooth  curves  with  a single  peak  for  Model  I. 

Curves  with  some  irregular  distributions  in  the  inner  portion  of  the 
wake  and  a single  smooth  peak  in  the  outer  portion  of  the  wake  occur 
for  Model  II.  For  Model  III  curves  with  double  peaks  and  a dip  were 
obtained.  The  height  of  the  peaks  decreases  in  the  downstream  direction 
except  at  near  stations  for  Model  I.  All  the  peaks  disappear  . . 
downstream,  and  the  rate  of  decay  for  Model  III  is  the  fastest  among 
the  three  models.  The  magnitude  of  the  axial  turbulence  intensity  in 
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the  wake  of  M>del  I is  relatively  lower  than  that  for  Model  III  which 
is  lower  than  that  for  Model  II.  The  radial  and  tangential  turbulence 
intensity  profiles  have  extra  ripples  after  the  first  peak  for  Model  I 
and  very  rough  distributions  for  Model  II.  For  Model  III  they  have  a 
higher  peak  after  the  first  peak  and  a dip  in  between  two  peaks  at 
near  stations.  The  maximum  value  of  these  quantities  decreases  down- 
stream for  the  three  models,  Model  III  has  the  highest  value,  and 
Model  I has  the  lowest  value  among  the  three  models.  The  radial  shear 
stress  profiles  are  smooth  curves  with  a single  peak  for  Model  I and 
curves  with  a peak  then  a big  dip  for  Model  III.  The  absolute  magni- 
tude for  Model  II  is  higher  than  that  of  Model  III.  Because  ;he  decay 
rate  of  these  quantities  for  Model  111  is  much  faster  than  for  the 
other  two  models,  the  magnitude  of  these  values  for  Model  III  are  even 
lower  than  for  the  drag  body  at  farther  stations  (X/D  = 20,  and  40). 

The  tangential  shear  stress  profiles  show  a high  peak  and  a big  dip  for 
Model  III.  The  magnitude  of  the  tangential  shear  stress  is  relatively 
low  for  Models  I and  II  as  compared  tr>  Model  III  which  exhibits  the 
highest  value  at  near  stations  but  it  also  has  the  highest  rate  of 
decay. 

Second,  for  the  cases  presented  here,  the  momentum  analysis 
suggests  that  the  turbulence  and  pressure  terms  are  of  non-negligible 
magnitude,  and  so  they  have  to  be  evaluated  n order  to  have  the 
drag  coefficient  exactly  zero.  This  i especially  necessary  for  the 
stations  where  the  turbulence  intensity  is  high  (Model  II)  and  the 
static  prepsure  variations  are  large  (Model  III).  According  to  t e 
results  presented  in  Section  "RESULTS  AND  DISCUSSION,"  Model  II  and 
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Model  III  have  very  reasonable  results,  although  the  pressure  and  tur- 
bulence terms  were  not  taken  into  account  in  evaluating  the  momentum 
analysis.  The  energy  analysis  suggests  that  the  resu  ts  presented  are 
good  in  general,  but  the  analysis  cannot  serve  as  a sensitive  detection 
of  the  accuracy  of  turbulence  measurements.  Nevertheless,  such 
analysis  gives  an  insight  into  the  overall  flow  structure  and  the 
mechanism  of  energy  transformation. 

Third,  the  wake  characteristics  of  the  blunt-tailed,  semi- 
streamlined  body  studied  by  Chevray  are  similar  to  the  wake  character- 
istics of  the  present  streamlined,  drag  body  (Model  I)  though  there 
remain  some  differences.  The  complex  vortex  system  near  the  tail  of  a 
drag  body  has  a significant  influence  on  the  turbulence  characteristics 
at  the  stations  near  the  body. 

Fourth,  the  wake  development  of  the  blunt-body  (disk)  which  was 
driven  by  a high-speed  central  jet  to  yield  a zero  momentum  wake  as 
studied  by  Naudascher  is  significantly  different  from  the  momentumless 
wake  development  of  a slender,  streamlined  body  driven  by  a peripheral 

jet . 

Fifth,  the  accuracy  of  the  measurements  in  this  work  is  within 
takeable  experimental  errors  according  to  the  check  of  the  momentum 
and  energy  balance  analyses. 

Finally,  since  there  are  few  available  results  regarding  moment un- 
less wakes  and  this  type  of  wake  provides  some  basis  for  t.ie  analysis 
of  submarine  detection,  further  development  and  study  in  this  area 
are  still  meaningful  and  practical.  Therefore,  It  is  recommended  that 
the  following  experiments  be  undertaken: 
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1)  Examine  the  effect  of  different  diameter  propellers  and  pro- 
pellers with  different  pitch  angles.  Also,  it  is  worthwhile  to  employ 
and  compare  a dual-propeller  case. 

2)  A swirling,  air  injection  model  should  be  developed. 
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APPENDIX  A 
SYMMETRY  OF  THE  WAKE 
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APPENDIX  A 


The  values  listed  in  the  Table  A-l  and  A-2  are  the  raw  axial 
turbulence  data  given  in  the  RMS  volt  units.  Since  the  RMS  value  is 
rather  insensitive  to  the  temperature  variations  (for  the  probes  and 
constant  temperature  system  applied  in  these  experiments) , it  can  be 
used  to  check  the  symmetry  of  the  wake  by  measuring  the  RMS  value  at 
points  which  are  at  the  same  distance  from  the  center  of  the  wake  but 
in  different  directions  (horizontal  left,  horizontal  right,  and 
vertical  downward) . 
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values  (In  volts)  recorded  by  straight  hot-wire  at  the  same  r but  different  directions:  (Model  I) 
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